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Abstract 

We report ab initio calculations of the electronic band structure and the phonon dispersion 
relations of the zincblende-type mercury chalcogenides (/3-HgS, HgSe, and HgTe). The latter have 
been used to evaluate the temperature dependence of the specific heat which has been compared 
with experimental data. The electronic band structure of these materials has been confirmed to 
have an inverted direct gap of the a-tin type, which makes HgSe and HgTe semimetallic. For (3- 
HgS, however, our calculations predict a negative spin-orbit splitting which restores semiconducting 
properties to the material in spite of the inverted gap. We have calculated the spin-orbit induced 
linear terms in k which appear at the Tg valence bands. We have also investigated the pressure 
dependence of the crystal structure and the phonons. 

PACS numbers: 63.20.-e, 63.20.dk, 63.20.D-, 68.35.bg, 65.40.Ba, 71.55.Gs, 71.70.Ej 
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I. INTRODUCTION 



The present work is part of an on-going effort to investigate, both theoretically and ex- 
perimentally, the temperature dependence of the specific heat ~ Cp of semiconductors 
and semimetals in the low temperature region in which strong deviations from the Debye 

law take place.- These deviations lead to a peak in Cp/T^ vs. T at a temperature T ^ 
0/20, where is the Debye temperature. Previous investigations have focused on diamond- 
like semiconductors and their zincblende and wurtzite analogs. -'^I'^i^'^ More recently the lead 
chalcogenides, with rock salt structure, have been investigated.- They, as well as the struc- 
turally related semimetals Bi and Sb, exhibit a remarkable dependence of their vibrational 
properties, and correspondingly C^, on spin-orbit (s-o) interaction.-i^ C„ calculations, and 
the corresponding Cp measurements, have also been reported for the related semiconductor 
(rhombohedral structure) GeTe, a material which grows with a large hole concentration that 
enables the investigation of the effect of such doping on the vibrational properties and on 
Cp (or Cv).— This work usually included measurements on samples prepared with differ- 
ent isotopic compositions, so as to investigate their effect on Cp , accompanied by similar 
calculations of Cy. 

Because of the calculated dependence of the aforementioned properties on s-o interaction 
for materials composed of heavy atoms (e.g. Bi, Pb) and structures related to rock salt, 
we decided to investigate similar materials that crystallize in the zincblende structure. The 
mercury chalcogenides appeared to be ideal for this purpose, although the toxic nature of 
mercury dissuaded us from investigating isotope effects. Moreover, the seven stable isotopes 
of Hg have low isotopic abundances and correspondingly prohibitive prices. These materials, 
however, have a number of hitherto not fully investigated electronic properties, related to the 
s-o interaction and to the gap inversion, which makes of them zero-gap semiconductors.— 
We have investigated some of these properties in the process of calculating ah initio the 
phonon dispersion relations and Cv 

HgS is usually found in nature as cinnabar, crystallized in the Se-like trigonal structure 
(a-HgS), a structure in which it usually also appears when grown in the laboratory. A 
small addition of Fe helps to crystallize HgS in the zincblende structure (/3-HgS), both in 
nature ((~1%) Fe), metacinnabar— ) and in the laboratory.— We use for the present work 
natural crystals of /3-HgS originating from Mount Diablo Mine, Clayton, CA (purchased 
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from John Betts, New York, NY 10025). A chemical microanalysis of two such crystals taken 
from different parts of the sample gave an Fe content of ~ 1.1 mass-% Fe. Additionally, 
we measured the specific heat of a commercially available polycrystalline fine powder of 
99.999% stated purity (Alfa Aesar GmbH & Co KG, Karlsruhe). HgSe and HgTe can also 
be found in nature (HgSe: Tiemannite, HgTe: Coloradoite) but they are rare. We used for 
the work reported here synthetic crystals grown by the Bridgman technique. Cinnabar is 
transparent in the red and has an optical gap of ~2.275 eV.— We shall not discuss it here 
any further since we confine ourselves to the chalcogenides with zincblende structure. It will 
be discussed in a later publication. The zincblende variety of HgS (/3-HgS, metacinnabar), 
which will be fully discussed here, has a black color and a zero or near zero gap. All crystals 
used here were found to be n-type by thermoelectric power measurements. 

The structure of this article is as follows: In Sec. [TTlwe discuss the procedures employed 
for the ab initio calculations of the electronic band structure, the phonon dispersion rela- 
tions and the heat capacity. Section IIIII presents some important results obtained in the 
process of performing the band structure calculations. We display the full LDA electronic 
band structure of /9-HgS, which has not been reported in the literature so far. We have 
found, however, three articles discussing the band structures of HgSe and HgTe which make 
it unnecessary to reproduce here our full calculated band structures for these materials. 
The oldest, by Bloom and Bergstresser— uses an empirical pseudopotential method which 
includes s-o interaction. The second one, by Overhof^, uses the KKR (Korringa-Kohn- 
Rostocker) method. In spite of not being self-consistent, these calculations already reveal 
the inverted, a-tin like band structure of these materials, as it does the self- consistent LMTO 
calculation for HgTe performed by Cade and Lee.— In the process of analyzing our band 
structure calculations, we realized that there is a considerable contribution of the core 5d 
levels of Hg to the s-o splitting at the F point of the valence bands. Like in the case of 
the copper monohalides, this contribution has a sign opposite to that of the splitting of the 
p valence electrons of the anion.— »^ For HgS, this 5d contribution overcompensates that 
of the 3p states of sulfur, thus reversing the corresponding valence states and changing the 
topology of the bands. This situation is similar to that found in CuCl, ZnO and CuGaS2.— 
The calculated s-o splittings for the three Hg compounds are given in Table [B This table 
also contains calculated values of the negative Fg - Fi gap and the lattice parameter 
obtained through energy minimization. Experimental data for these parameters are also 
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given. 

Another interesting aspect of the band structures of the HgS, HgSe and HgTe is the 
appearance of hnear terms in the wavevector k at the Fg point of the valence bands. Such 
terms, of relativistic origin, were predicted by G. Dresselhaus for zincblende-type materials 
on symmetry grounds.— They arise from s-o coupling plus the lack of inversion symmetry 
and involve k-p interactions between the valence bands and the core d electrons.— Very few 
ab initio calculations of these terms, which are important in present day spintronics, have 
been performed and the Hg chalcogenides should be ideal materials for such calculations. 
We thus present enlarged versions of the valence band structures around F and use them to 
evaluate the coefficients of these linear terms. The results are presented in Table [ITl 

Section HVl is devoted to the ab initio lattice dynamical calculations and comparison with 
inelastic neutron scattering and Raman data. For the purpose of interpreting the specific 
heat results, we also present the corresponding one-phonon densities of states and their de- 
composition in Hg and anion contributions. The possible contributions to the linewidths 
of the Raman phonons, anharmonicity, electron phonon interaction and isotopic disorder is 
also discussed. For completeness, we also present the density of two-phonon states (sum and 
difference) corresponding to a total zero k-vector. These results are useful for the interpre- 
tation of two-phonon Raman scattering data; they are compared with spectra available in 
the literature for HgSe. 

Because of the lack of theoretical information concerning the phase transitions of the 
Hg chalcogenides under pressure, we have also performed calculations of the enthalpy of 
HgS and HgTe versus pressure for the zincblende, rocksalt and cinnabar phases. From the 
calculated dependence on pressure p, the values of p at which transitions occur are obtained 
and compared with experimental results. As a by-product we report calculated values of the 
Griineisen parameters of the LO and TO phonons of HgTe at the F-point of the Brillouin 
zone. This work is presented in Sect. |Vl 

Section IVII discusses the specific heats and their temperature dependence around the 
maximum of C^/T^- The ab initio calculations, which have been performed with and without 
s-o interaction and using either the LDA or the GGA approximation, describe rather well 
the experimental results. The inclusion of s-o interaction lowers the maximum of Cy/T^ 
by about 8% in the three compounds with zincblende structure. Since this effect has a 
sign opposite to that found for the rock salt structure semiconductors (e.g. PbS), we tried 
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to ascertain whether the sign of the effect depends on crystal structure. For this purpose, 
we calculated the lattice dynamics and of HgS in the pressure region where the rock 
salt phase is the stable one. We found a very small effect of the s-o interaction at the 
maximum of C^/T'^, corresponding to a decrease by less than 1%. Section rVIII summarizes 
the conclusions. 

II. THEORETICAL DETAILS 

In this paper we study the vibrational and thermal properties of HgX (X=S, Se and Te) 
by three different approaches: an experimental one and two different calculations based on 
different implementation of density functional theory.— For the theoretical part, we have 
used two different computer codes: ABINIT and VASP.— In order to describe the ex- 
change and correlation effects we have used two well known approximations, LDA (local 
density approximation)^'^ and PBE-GGA (Perdew, Burke, Ernzerhof generalized gradient 
approximation) .— 

In the ABINIT calculations, the wavefunction is expanded in plane waves and a norm- 
conserving pseudopotential approximation to the free ion potential is employed, with the 
Teter-Pade parametrization for the exchange and correlation.— In particular, we did use 
the Hartwigsen-Goedecker-Hutter HGH pseudopotentials, where scalar-relativistic and SO 
terms are combined, such that they fully represent nonlocal pseudopotentials of separable 
form.— 1^1^ A 60 Ha kinetic energy cutoff was applied and a Monkhorst-Pack grid of 6x6x6 
was used in order to describe any electronic quantity within the Brillouin zone. Electronic 
and vibrational properties were checked with respect to these parameters to ensure conver- 
gence. The lattice dynamics was calculated by using an implementation of density functional 
perturbation theory.— Vibrational frequencies and the vibrational density of states were cal- 
culated after the diagonalization of the interatomic force constants (IFC). These were used 
to obtain the thermal properties within the harmonic approximation. An original grid of 
12x12x12 within the Brillouin zone, with a total of 72 q points, was used to calculate the 
corresponding set of IFCs. Afterwards, a finer grid was employed to determine the thermal 
properties, obtained by performing a Fourier fit to the calculated IFCs (for example in order 
to calculate C^, a set of 32768 q points was used, which was reduced by symmetries to 897 
points). 
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The Vienna ab initio simulation package (VASP) is a density functional theory imple- 
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mentation with pseudopotentials used to obtain the electron wave functions [see Ref. 
and references therein]. The exchange-correlation energy was taken in the local density 
approximation (LDA) with the Ceperley Alder prescription. Calculations with the general- 
ized gradient approximation (PBE-GGA) were also performed.— The projector-augmented 
wave (PAW) scheme was adopted and the semicore 5d and 6s electrons of Hg were included 
explicitly in the calculations.— The set of plane waves used extended up to a kinetic energy 
cutoff of 370 eV for HgS and 310 eV for HgTe and HgSe. This large cutoff was required 
in order to achieve highly converged results within the PAW scheme. We used a dense 
Monkhorst-Pack grid for BZ integrations in order to ensure highly converged results in the 
three compounds (to about 1-2 meV per formula unit). We also used an accurate algorithm 
during the calculations which allowed us to obtain highly converged forces for the calculation 
of the dynamical matrix. At each selected volume, the crystal structure was fully relaxed to 
its equilibrium configuration through the calculation of the forces on atoms and the stress 
tensor.— In the relaxed equilibrium configuration, the forces are less than 0.002 eV/A and 
the deviation of the stress tensor from a diagonal hydrostatic form is less than 0.1 GPa. 

Highly converged results on forces are required for the calculation of the dynamical matrix 
using the direct force constant approach (or supercell method).— The construction of the 
dynamical matrix at the P-point of the BZ is particularly simple and involves separate 
calculations of the forces in which a fixed displacement from the equilibrium configuration 
of the atoms within the primitive unit cell is considered. Symmetry helps to decrease the 
number of such independent distortions, thus reducing the computational effort in the study 
of the analyzed structures considered in this work. Diagonalization of the dynamical matrix 
provides both the frequencies of the normal modes and their polarization vectors which 
allow us to identify the irreducible representation and the atomic character of the phonons 
modes at the relevant k point. Supercell calculations were carried out in order to obtain 
phonon dispersion curves. Starting with a 2x2x2 supercell the direct method delivers 
correct harmonic phonon frequencies at the special wave vectors of the Brillouin zone which 
are compatible within the selected supercell size. If the interaction range would cease at 
the selected supercell size, the corresponding phonon frequencies would provide a good basis 
for the interpolation to all wave vectors. The LO/TO splitting cannot be automatically 
included when using the direct method, but the dynamical matrix can be supplemented 
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with a non-analytical term which depends on the Born effective charge tensors and the 
electronic dielectric constant.— The thermodynamic properties of the crystals are, to a large 
extent, determined by phonon dispersion and can be satisfactorily described in most of the 
cases within the quasiharmonic approximation, as discussed in Ref. |36. 



III. ELECTRONIC BAND STRUCTURE 



Figure 1 displays the band structure of /3-HgS as calculated with the ABINIT code using 
the LDA approximation (see Sect. [IT]) with inclusion of s-o interaction. The band structure 
shown in this figure was calculated with the energy optimized lattice parameter ao = 5.80 
A. The s-o splitting of the Fis valence bands (-0.18 eV) cannot be resolved in this figure 
because of the wide energy scale it covers. It will be discussed later in connection with Fig. 
2. The wide energy scale, however, allows us to observe the position of the lowest, sulfur 
3s-like Fi valence band, which is not affected by s-o interaction, at about -12 eV. Above 
this band, the 5d bands of Hg appear. Without s-o interaction they split into an upper 
quadruplet (~-6eV, F12) and a lower sextuplet (F15). s-o interaction does not split the F12 
states but it splits the F15 sextuplet into a doublet F7 and a quadruplet Fg , the former being 
above the latter. The splitting shown in Fig. 1 is 2.1 eV , the same as the corresponding 
atomic splitting of Hg.— 




FIG. 1: Electronic band structure of /3-HgS calculated with the ABINIT code in the LDA approx- 
imation (see SectUT]) and with s-o interaction using the lattice parameter oq = 5.80 A obtained by 
total energy minimization. 
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At -0.58 eV we find an 5-like Fi state which is composed mainly of 6s states of Hg. 
This state corresponds to the lowest conduction band for most tetrahedral semiconductors 
(exception: a- tin^i and, as discussed below, HgSe and HgTe) above the Fis valence bands. 
In order to resolve the details of the band structure around the F point we replot in Fig. 
|2] the data from Fig. [T] with an expanded energy scale. This figure shows the inverted 
s-o interaction at F: the F7 doublet is above the Fg quadruplet. This fact arises from 
the contribution of the negative s-o splitting of the F15 states of Hg (Fig. 1) which 
overcompensates the positive splitting of the S 3p contribution. This effect is similar to that 



observed for CuCl (Refs. [l8|,19|) but for Hg an additional complication arises, related to the 
a-tin-like structure. Counting bands we see that the Fg bands should be occupied whereas 
the s-o split F7 should be empty. Hence, according to Figure 2, undoped (intrinsic) /3-HgS 
should have an indirect energy gap of ~0.15 eV and a direct one of 0.18 eV . This is in 
reasonable agreement with the gap observed by Zallen and Slade (0.25 eV), although these 
authors did not realize that the s-o splitting of /3-HgS is negative.— 




Wavevector 



FIG. 2: Electronic band structure of /9-HgS calculated with the ABINIT code in the LDA approxi- 
mation (see Sect. HI]) using the lattice parameter ao = 5.80 A obtained by total energy minimization 
as in Fig. [T]but with the energy scale enlarged so as to display the details of the top valence bands 
and the lowest conduction band around F. The Fg - r7 separation corresponds to the negative s-o 
splitting Aq. Fg - Fg is the Eq gap. 

In Table H] we have listed the negative gaps Eq between the valence bands and the Fg Hg 
5s-like conduction band calculated and observed for the three Hg chalcogenides. We notice 
in the table that the calculated positions of the Fg bands are lower then the experimental 
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TABLE I: Values at k=0 of the spin- orbit splitting Aq, the negative p-s gap Eq (in eV) and the 
lattice constant cq (in A) as calculated with various codes and experimentally determined. Notice 
the negative value of Aq calculated by our three methods for HgS. It is due to an admixture of p- 
wavefunctions of the cation with 5d of Hg (the latter have a negative contribution similar to results 
found for CuCl, ZnO, CuGaS2)- The differences of Aq for the different computational methods are 
due to minor differences in the compensating p and d contributions. The negative Aq of /3-HgS 
has already been mentioned recently.— 



HgS HgSe HgTe 





VASP-LDA 


-0.111 


-1-0.23 


-1-0.783 




A TITIVTT'T"' T T "\ A 

ABiJ\li-LDA 


-0.18 


-1-0.34 


-1-0.81 


Ao (eV) 


VAbP-GGA 


-0.091 


-1-0.25 


-1-0.751 




(Jverhof 




-1-0.48 


-1-1.13 




experimental 


-0.25 (Ref. 12) 


+0.38 (Ref. 41) 


0.9 (Ref. 41) 




VASP-LDA 


-0.573 


-1.18 


-1.025 




ABINIT-LDA 


-0.58 


-1.26 


-1.15 


^0 (eV) 


VASP-GGA 


-0.483 


-1.07 


-1.113 




Overhof 




-0.4 


-0.3 




experimental 


-0.15 (Ref. 12) 


~-0.4 (Ref. 41) 


-0.3 (Ref. 41) 




VASP-LDA 


5.825 


6.065 


6.433 




ABINIT-LDA 


5.80 


6.054 


6.442 


ao (A) 


VASP-GGA 


5.999 


6.255 


6.633 




ABINIT-GGA 


5.982 


6.250 


6.656 




Overhof 




6.075 


6.45 




experimental 


5.841(~295 K) (Ref. 42) 


6.085 (293K) (Ref. 41) 


6.453 (Ref. 41) 



ones by 0.43 eV for HgS and 0.85 eV for HgSe and HgTe. It is natural to attribute this 
fact to the so called "gap problem" inherent to the LDA calculations.— Recent GW show 
that the Fi state of the HgX materials should be raised by about 0.75 eV.— If we add 
this energy to the values of Eq displayed in Table [H the Eq gaps of HgSe and HgTe remain 
negative (~ -0.35 eV) whereas that of /3-HgS becomes positive (~ -1-0.4 eV). Thus /3-HgS 
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would be a conventional semiconductor except for the negative Aq. As already mentioned, 
we do not reproduce here the full band structures of HgSe and HgTe because they have 
already been presented in the literature.— >i^>i^ We have given a few important parameters 
in Table [U also for HgSe and HgTe. An interesting fact is their positive s-o splitting Aq 
as opposed to the negative one found for HgS. The situation is similar to that found for 
the copper halides: CuCl has a negative splitting whereas for CuBr and Cul it is positive.— 
This fact has been attributed by Shindo and Kamimura to the negative contribution of the 
copper 3d electron to the valence band splitting.— As already suggested above, a similar 
explanation holds for the Hg chalcogenides: the negative splitting of the 5d electrons of Hg 
contributes, through hybridization, an amount to the valence splitting sufficient to reverse 
the 3p splitting of S but neither the 4p splitting of Se nor the Sc/ of tellurium. According to 
Table I, agreement between experimental and calculated values of Aq is good. Notice that, 
because of the positive s-o splitting, HgSe and HgTe are zero gap semiconductors, like a-tin. 

In the process of performing these calculations we also took a look at the terms linear in 
k which appear in the valence bands around F as a result of s-o interaction coupled with the 
lack of inversion symmetry. The existence of these terms was first reported by Dresselhaus^ 
and numerical calculations were performed by Cardona, Christensen and Fasol^, by Surh, 
Li and Louie^ and more recently by Kim et al.—. These terms arise in perturbation theory 
from s-o interaction between the valence bands and the d core electrons (ri2 states) and back 



to the valence band via k-p.— In Ref. |22| ab initio calculations of these terms were presented 
for many III-V and II- VI zincblende-type semiconductors but not for the Hg-chalcogenides. 
These calculations may have been vitiated by "gap problems" but concerning the linear 
terms in k any such errors should be small because they do not involve the p-s gap. In Ref. 



22 



a semiempirical expression was given for the coefficient Ck of these linear terms. For the 



Hg-chalcogenides compounds it can be written as (see Eq.(7.6) of Ref. 
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Where A^ is the d s-o splitting of the cation (for Hg 2.2 eV) and the energy denominator 
represents the gap between the Fg valence band and the outermost d-states of the cation. 
In EqlU we have neglected the contributions of the anion d electrons since they are absent 
for HgS and should be small for HgSe and HgTe. Equation [1] allows us to estimate Ck = 
-128 meVA for the Hg-chalcogenides. 
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We now give some details concerning the definition of Ck- Tliis parameter lias been 
chosen such that the splitting for k along [100] is ±Cfc k (the convention used for the sign is 



that of Ref. 



221 ). each one of the two split states being a doublet for this k direction. For k 



along [111] two of the four Fg states are not affected by linear k terms, whereas the other two 
split like ±V2Ck A; . For k along [110] the Fg states split into 4, two like ±Ck k{l + 73/2)^/^ 
and the other two like ±Ck k{l - v^/2)V2_.2i^ 

We display in Fig. [3] the valence bands of HgSe calculated in a small energy and k range 
so as to be able to observe the linear terms and thus determine C^- The calculations were 
performed for k along [100] and [111] using the ABINIT code in the LDA approximation. 
Values of Ck were obtained through fits with a polynomial containing a linear and a quadratic 
(effective mass) term. Notice the splitting pattern mentioned above: for k along [100] the Fg 
quadruplet splits into two doublets. For k along [111] a doublet does not split and exhibits 
no linear term (notice the small value of the linear fitting coefficient). The other doublet 
splits into ±-\/2Cfc k. Similar patterns are found for HgS and HgTe. We shall not present 
here the corresponding figures. Instead we shall list all of the values of Ck determined for 
our three materials with the ABINIT and the VASP codes. 




FIG. 3: (colour online) Fg valence bands of HgSe showing the linear terms in k along the [100] and 
the [111] directions, as calculated with the ABINIT-LDA code. The colored curves represent fits to 
a quadratic polynomial, the black ones connect calculated points as a guide to the eye. The values 
of the fit coefficients are given in the inset. The zero of energy lies at the valence band maximum 
at the Fermi energy. 

We display in Table [TTl the values of Ck for the three materials as found with the two 
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TABLE II: Values of the coefficients of the Unear terms around the Fg states (in meVA ) of the 
valence band structures of /3-HgS, HgSe and HgTe as calculated with the VASP and ABINIT codes 
through a quadratic fit of the eigenvalues (see Fig. 3). 





/3-HgS 




;Se 




;Te 




Ck [100] 


Ck [111] 


Ck [100] 


Ck [111] 


Ck [100] 


Ck [111] 


VASP- GGA 


-147 


-140 


-145 


-146 


-145 


-137 


ABINIT-LDA 


-144 


-174 


-160 


-170 


-148 


-136 



computational codes from the 



as that of Ref. 
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100] and the [111] splittings. The sign convention is the same 
and Ref. |46|. These values are rather close to that estimated with Eq. [1] 
using for the p-d gap 6 eV and 2.2 eV for the s-o splittings of the electrons of Hg {Ck = 
-128 meVA). We have not found experimental values of Ck for these materials. 



IV. LATTICE DYNAMICS 



A few Inelastic Neutron Scattering (INS) determinations of the phonon dispersion rela- 
tions of /J-ngS^"^, HgSe^ and HgTe^ have been published. They have been compared with 
semiempirical calculations which use fitted force constants. Here we compare the measured 
points with ab initio calculations (ABINIT or VASP) and examine the effect of s-o interac- 
tion on these calculations. We start with /3-HgS and display in Fig. Hlthe phonon dispersion 
relations calculated with s-o interaction (solid line) and also without (dashed line). The 
effect of the s-o interaction (see inset for the TA phonons) is rather small, however it is 
possible to see that for the TA phonons, those mainly responsible for the C^ to be discussed 
later, the frequencies without s-o are about 5% lower than with s-o. 

An interesting feature of Fig. H] is the strong upwards bending of the optical bands away 
from k = 0, the opposite of what happens for most tetrahedral semiconductors. This effect 
is due to the large mass difference of Hg and S: at the edge of the zone the "effective" 
mass which determines the frequency is basically that of sulfur, whereas at k = it is the 
somewhat smaller reduced mass. In order to obtain the calculated up-bending one must 
also invoke a strong force constant connecting the anions. The upward bending is reflected 
in the corresponding density of phonon states, shown in Fig. [5], together with those of HgSe 
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Wavevector 



FIG. 4: (color online) Phonon dispersion relations of /3-HgS calculated with the ABINIT-LDA code 
with (solid line) and without (dashed line) s-o coupling. The circles were obtained with INS^. The 
LO and TO phonons at k (red) squares were obtained by Raman scattering.— The frequency 
scales in the insets have been expanded to show the effect of the s-o around X. 

and HgTe to be discussed later. A direct consequence is that the TO Raman modes do 
not overlap the DOS continuum while the LO mode does. The width of the former should 
therefore not be affected by elastic scattering due to defects, such as isotopic fluctuations, 
whereas the LO-modes, overlapping the background, should. This situation is opposite to 



that discussed in Ref. 



51 



for ZnSe, where isotopic disorder is shown to broaden the TO 
Raman phonons but not the LO phonons. 

The up-bending of the dispersion relations are also important when considering confine- 
ment in quantum dots. These effects have been discussed for the case of CdS, a material for 
which the the TO bulk bands bend up but the LO bands bend down.— No such experiments 
are available, to the best of our knowledge, for /5-HgS but quantum confinement in dots of 
this material should produce an increase in frequency as opposed to most other related ma- 
terials. Up-bending is also found in the dispersion relations of HgSe and HgTe (see below), 
but not as strong as in the case of (3-B.gS because the difference between anion and cation 
masses is not as large. A canonical example of up-bending phonon dispersion relations is 
found in rock salt structure PbS. In this material the LO-phonons are found at 205 cm~^ 
at the r point of the BZ. A direct consequence of the up-bending is the fact that the cor- 
responding peak in the second order Raman spectrum is not centered at 2x105= 210 cm~^ 
but somewhat higher (228 cm~^ at 12K).— In contrast, for GaAs, a typical down-bending 

13 



case, the two LO-phonon peak is shifted down by as much as 5 cm~^ with respect to twice 
the frequency of the LO phonons at F.— Effects of the upwards bending on the phonon hne 
widths will be discussed below. 




FIG. 5: (Colour online) Phonon density of states of /3-HgS, HgSe and HgTe, (a), (b), and (c), 
respectively, as obtained with the ABINIT code used for Fig. HI We have plotted the total DOSs 
and also their projection on the Hg and the anion atoms. Note that in the case of /3-HgS the 
acoustic phonons correspond to nearly pure vibrations of Hg, the optic ones to nearly pure S 
vibrations. The (blue) vertical arrows mark the positions of the optical phonons at F. 

Figures El and [7] display the phonon dispersion relations of HgSe and HgTe, respectively, 
as calculated with the ABINIT-LDA code. In Fig. [6] we only show the acoustic phonon 
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branches, and the corresponding INS data because ab initio calculations, using the VASP- 
ASA code, have already been reported by Lazewski and coworkers.— The main purpose of 
this figure is to visualize the effect of the s-o interaction which slightly raises the frequencies 
of the TA modes and, as we shall see below, lowers the calculated maximum of Cr^/T^ . 



The optical branches calculated in Ref. |55| for HgSe also bend upwards when moving 



away from F, but not as much as in Fig. HI This is expected because for HgSe anion and 
cation masses differ less than in the case of /3-HgS. Correspondingly, and as shown in Fig. [5l 
the optical phonon bands of HgSe are narrower than those of /3-HgS. Moreover, the projected 
densities of states of these phonons correspond about 90% to Se and 10% to Hg. 




FIG. 6: ah initio calculations of the acoustic dispersion relations of HgSe performed with the 
ABINIT-LDA code, with (solid line) and without s-o (dashed line) interaction. The circles repre- 



sent INS data taken from Ref. l55|. We also performed calculations with the VASP-GGA code. The 
LA modes were found to be slightly lower in frequency (82 cm^^ at L, 91 cm~^ at X) as expected 
from the larger oq (see Table Hj). The TA modes barely changed. 

Figure [7] displays the full acoustic and optic dispersion relations for HgTe together with 
INS data. In this case the upward bending of the optical bands is rather small, as corresponds 
to cation and anion masses which differ by less than a factor of two. Because of the weak 
upward bending, the optical DOS band is even narrower than for HgSe (see Fig. [5]). 

For the sake of completeness we include a figure with the two-phonon DOS which applies 
to second order optical processes in which either two phonons are created (sum processes) 
or one is created and the other destroyed (difference processes). After multiplication by the 
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Wavevector 



FIG. 7: Phonon dispersion relations calculated with the ABINIT-LDA code for HgTe. The circles 
represent data from INS measurements Raman spectroscopy gives TO(r) = 116 cm~^ and and 
LO(r) = 139 cm-^^ 

appropriate matrix elements (which we have not calculated in this work) the two phonon 
DOS should describe the corresponding IR and Raman spectra. Even without multiplication 
by the matrix elements, structure in the two-phonon DOS should help to interpret the 
corresponding spectra. 

There is only scarce information about the second order Raman and IR spectra of the 
materials under consideration. We have already mentioned the 2L0 structure seen for /3- 
HgS, which is slightly upshifted with respect to twice the frequency of the LO (F) phonons. 
Szuszkiewicz et al. have reported Raman spectra for HgSe with structures which can be 
attributed to two phonon sums.— A rather detailed assignment is now possible using the 
DOS of Fig. [HI The resolution of their measurements is better than 2 cm~^. The FWHM 
of the observed TO phonon is ~5cm~^. After correcting for instrumental resolution^, the 
FWHM becomes ~4.2 cm~^. We shall discuss below a possible isotope disorder contribution 
to this width.— The temperature at which the spectra were measured (90K) most likely 
precludes the observation of difference processes. The measurements were performed in 
back-scattering using a [110] surface of cleavage. Under these conditions, the selection rules 
prevent the observation of single LO phonons. The single TO mode is seen at 133 cm~^. We 
show in Table UTTl the frequencies of the observed "two phonon" peaks and their assignment 
to structures in Fig. [HI We also have made an attempt, as given in Table IIII| to assign 
these structures to the dispersion relations calculated for HgSe. The exact overtones of 
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FIG. 8: (color online) Sum ((black) solid line) and difference ((red) dashed line) of the two-phonon 
DOS of /3-HgS, HgSe and HgTe. The wavevectors of the two phonons are taken to be equal (but 
opposite), as required by optical spectroscopies. 

the phonons at T [2T0(r)= 264 cni-^ and 2L0(r)=330 cm ^ ] are not clearly observed: 
2T0(r) corresponds to a threshold and not a peak (see Fig. [8]), 2L0(F) should be slightly 
shifted upwards because of the up-bending and may overlap with the broad band that we 
have assigned to 2L0(X), as indicated in Table UTTl 
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TABLE III: Structures observed in the second order Raman spectrum of HgSe (Ref. 1571 ) compared 
with structures seen in the two-phonon DOS calculated on the basis of the band structure partly 
depicted (only the acoustic bands) in Fig. [71 



experiment (cm ^) calculated (cm ^) assignment 

43 50, 65 2 TA 

110 110 LA + TA 

120 120 LA + TA 

125 130 LA + TA 

160 (weak shoulder) 

195 199, 200, 210, 220 O + TA 

260 260, 245, 275, 280 2 TA, LO + TA 

350 230 (exactly 2 LO(r)) 2 LO close to F, 2 LO(X) 

370 370 2L0 (L,S) 



A. Line widths of F-phonons and isotope disorder contribution to their broaden- 
ing 

We have found reliable data for these linewidths in the cases of HgSe^ and HgTe^. The 

for HgTe is 1 cm~^ (FWHM). Correspondingly, no resolution 
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resolution reported in Ref. 
correction is needed for the measured linewidths which are 5 cm~^ for the TO (F) phonon 
and 4.5 cm~^ for LO (F), at 1.6 K.— In the work of Szuszkiewicz on HgSe the resolution 
is quoted to have been better than 2 cm~^, standard deconvolution:^ of the observed line 
width of the TO phonons (5 cm~^) brings the measured FWHM to about 4.2 cm~^. These 
materials being semimetallic, there should be some contribution of transitions from the 
occupied to the empty valence bands to the phonon linewidths. An investigation of this 
contribution has not been performed. The temperature dependence due to anharmonic decay 
into two phonons^, has not been investigated for these materials, neither experimentally nor 
theoretically. Nevertheless extant data, even those obtained at the lowest temperatures, are 
affected by anharmonic effects related to the zero-point vibrations.— These zero point effects, 
however, are typically less than 1 cm~^ (FWHM) for most zincblende-type semiconductors.— 
We thus surmise that the electron phonon coupling across the Fg valence bands by either 
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TABLE IV: Isotopic fluctuation parameter g (x 10^) for the elements under consideration, either 
as found in nature or prepared as equal mixtures of the two isotopes whose masses are given in 
brackets. 



isoptope composition 


Hg 


S 


Se 


Te 


natural 


6.2 


9.2 


7.9 


3.8 


50-50% 


(198-204) 22 


(32-34) 92 


(76-80) 67 


(124-130) 56 



the LO or the TO F-phonons must give a sizeable contribution to their FWHM. Before we 
assign to electron-phonon coupUng the difference between the observed FWHM (~5 cm"^) 
and the expected anharmonic effect (~1 cm~^) we must, however, estimate the contribution 
of the isotopic disorder in materials grown from the natural elements. This can be done 
with the expression:— 

= ^^[gngN^^^uo) leoHgP + fi-ch iVdChlt^o) leochH (2) 

where Fj represents the contribution of isotopic disorder to the FWHM of either the TO 
or the LO phonon, (^Hg and (^ch the isotopic mass fluctuation parameter of either Hg or the 
chalcogen, respectively: 

_ (^Hg,Ch) - (MHg,Ch)^ 

^Hg.Ch - 7TT l-^j 

WwHg,Ch/ 

In Eq. [2] A^^Hg and Ndch are the densities of one-phonon states projected on Hg and the 
chalcogen, respectively, while eoHg and eoch are the corresponding eigenvector amplitudes 
for the phonons at F. We list in Table IIVI the values of g for the natural elements under 
consideration and for artificial isotopic mixtures conceived so as to maximize g taking into 
account the price of the isotopes. 

We thus conclude that the observed FWHM of TO(F) and LO(F) phonons in HgSe and 
HgTe, between 3 and 4 cm~^ after subtracting anharmonic effects, must be due to electron- 
phonon interaction between full and empty valence bands. The effect can be compared to 
that observed in heavily doped Si and Ge. For Si with 6x10^^ cm~^ holes a contribution of 
about 5 cm-i to the FWHM is found.-^ 
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TABLE V: Parameters needed to evaluate the disorder induced broadening of LO (F) phonons 
with Eq. [2] and the isotope disorder induced broadenings obtained by this procedure. 





NdRg 


Ndch 


eoHg 


eoch 


Ti (natural) (cm ^) 


7i (50-50%) (cm-i) 


/3-HgS 





0.015 


0.15 


0.99 


0.02 


0.2 


HgSe 


0.03 


0.085 


0.40 


0.93 


0.06 


0.5 


HgTe 


0.02 


0.05 


0.53 


0.85 


0.02 


0.2 



V. PRESSURE EFFECTS ON THE STRUCTURE AND LATTICE DYNAMICS 



This investigation was triggered by our interest in calculating the specific heats of the 
Hg chalcogenides crystallized in the rock salt structure (see next section). These materials 
undergo phase transitions at relatively low pressures.— While there is considerable exper- 



imental data on these transitions (see Ref. 



62l ). no ab initio calculations seem to have 



been reported. We have performed enthalpy minimizations vs. pressure for HgS and HgTe 
considering the three most conspicuous structures of these materials: cinnabar, zincblende 
and rock salt. Orthorhombic and bcc structures have also been observed experimentally for 
HgTe but we have not included them in our calculations. 

We show in Fig. [9]the pressure dependence of the enthalpy of the three phases mentioned, 
as calculated for HgS. The cinnabar phase (a-HgS) is found, within estimated error, to have 
the same enthalpy as the zincblende phase (/3-HgS), a fact which may explain the ease to 
obtain both phases in the laboratory, although the growth of the zincblende phase requires 
the addition of a small amount of Fe.— 

Notice, however, that the zincblende phase is calculated to be the stable one at zero 
pressure, but with an enthalpy only about ~0.0003 Ha below the cinnabar phase. We 
do not believe that this minute enthalpy difference is significant. The cinnabar phase is 
calculated to be stable between 0.7 GPa and 22.2 GPa, at which pressure the NaCl phase 
should appear. This is in reasonable agreement with reports of a phase transition at 20.5 
GPa.— In the next section we shall present calculations of the specific heat of the rock salt 
phase of HgS in the stable pressure region, so as to ascertain the effect of s-o interaction on 
this phase. 

In Figure [TO] we present similar calculations for the enthalpy of HgTe vs. pressure. 
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FIG. 9: (color online) Enthalpy of the three phases of HgS calculated with the ABINIT-LDA code 
with s-o interaction. Cinnabar and rock salt phases are basically degenerate (within any reasonable 
error estimates) at p=0 although formally their enthalpies cross at ~ 0.7 GPa. The cinnabar — > 
rock salt transition is predicted to occur at ~22 GPa. Experimental data place it at 20.5 GPa.— 



0.010 




-0.010 



-2 2 4 6 
Pressure (GPa) 



FIG. 10: Enthalpy of HgTe calculated with the ABINIT-LDA code including s-o splitting. The 
arrow signals the calculated zincblende — > cinnabar transition which is found experimentally to 
take place at 1.5 GPa.— 

The calculations were also performed without s-o interaction and we found that the 
presence of this interaction slightly improves agreement with experimental results which are 
1.5 GPa for the zincblende cinnabar transition (as compared with 1.3 GPa in Fig. [TOl 
1.0 without s-o interaction. 8 GPa for the cinnabar rock salt transition has been also 
observed.— For reasons unbeknownst to us it is not reproduced in our calculations. 
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As a complement to this pressure investigation we calculated the pressure dependence of 
the mode frequencies of the LO(r) and TO(r) phonons for zincblende HgTe. The results, as 
obtained with the ABINIT-LDA code with and without s-o interaction, are shown in Fig. 

m 




FIG. 11: Pressure dependence of the frequencies of the LO(r) (upper curves) and TO(r) phonons 
(lower curves) of HgTe calculated with the ABINIT code with and without s-o interaction as 
indicated. 

They can be expressed in terms of the usual Griineisen parameters using as a bulk mod- 
ulus B= 42.2 GPa:ii 

dlnuLo din uto 

= -^I^aT = 1-0' = -^k^ = '-^^ 

using Fig{TT] it is also possible to determine a Griineisen parameter for the transverse 

effective charge e* as defined with the equation: 

2 2 _47riVe*2 

'^LO - '^TO - (5) 

The calculation of a Griineisen parameter (i.e. a logarithmic derivative with respect to 
the volume), like in Eq. |4]) for e* from the data of Fig. [11] requires, beside the bulk modulus, 
the knowledge of the Griineisen parameter of the electronic dielectric constant eoo-— 
For lack of better information we take a value of = +1 , typical of that found for most 
tetrahedral semiconductors^ and find for the Griineisen parameter of e*: 7*= -0.65. This 
negative sign implies that e* increases with increasing volume. It is usually interpreted 



22 



to signal an increase in ionicity with increasing volume, a result that is common to most 
zincblende-type semiconductors with the exception of SiC (see Table IV in Ref. l65l ). 



VI. SPECIFIC HEAT 

The specific heat Cp of the crystalline samples was measured for the three chalcogenides 
between 2 and 280 K with a PPMS system (Quantum Design, 6325 Lusk Boulevard, San 
Diego CA (Quantum Design, 6325 Lusk Boulevard, San Diego, CA), as described before.-*^ 
The powder samples of /?-HgS were sealed into glass fiasks under 900 mbar ^He gas to 
ensure thermal contact at low temperatures. The heat capacities of such samples were 
measured using a quasi-adiabatic heat pulse method with a home-built Nernst calorimeter 
similar to that described Schnelle and Gmelin.~ The heat capacities of the glass sample 
container, the sapphire platform and a minute amount of Apiezon N vacuum grease used to 
thermally couple the sample to the platform were measured separately and subtracted. 

The specific heat of HgSe above 80 K has been measured by Gul'tyaev and Petrov-^. Low 
temperature specific heats of HgTe have been determined by Collins et al. and subsequently 
for /3-HgS, HgSe and HgTe by Khattak and coworkers.— Our data for HgSe and HgTe 
agree very well with those measured previously. We find, however, very large differences 
at temperatures around 7.5 K where the maximum in Cp/T^ occurs for our /?-HgS data 
compared to those of Khattak et al., our data being a factor ~2 larger. 

Calculations of were performed by integrating the calculated phonon DOS with the 
standard expression (we have already mentioned that in the region of interest Cp ~ C^). 
The calculations were performed both, with the ABINIT-LDA and the VASP-GGA codes, 
with and without s-o interaction. The effect of s-o interaction was basically the same for 
the two codes: the s-o interaction lowered the maximum of C^/T^ by about 10%. We shall 
present here the results obtained with only one of the codes (except for /3-HgS) but both, 
with and without the s-o interaction. 

We show in Fig. WMa) the experimental values of C^jT'^ for two sample of /3-HgS, a 
natural crystal and a high-purity powder. The natural crystal was found by chemical micro- 
analysis to have an Fe content of 1.1 mass-%. For comparison we also display measurements 
by Khatak and coworkers.— The upturn at low temperatures observed for the natural crys- 
tal can be ascribed to a linear term in the specific heat (Sommerfeld term) due to metallic 
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behavior. To correct for this term we subtracted a hnear contribution 7=20 mJ/molK2) 
from the experimental data. Corrections for an FeS content using the data by Kobayashi et 
ali^ lead to a negligible change of the data and is not shown here. 

Fig. TMh) displays calculated VASP and ABINIT values of CJT'^ for /3-HgS obtained 
with and without s-o interaction. The agreement between theoretical and experimental data 
is reasonably good. 




FIG. 12: (color online) (a) C^/T^ measured on a mineral and a powder sample of /?-HgS. Data by 
Khattak et al. are also given for comparison.— The solid line represents the mineral data minus a 
linear contribution as indicated in the inset, (b) Same mineral data as before compared with the 
results of four calculations as indicated in the inset. 

Figure [13] shows Cp/T^ measured by us and the values calculated from data by Gul'tyaev 
and Petrov^, compared with ABINIT- LDA calculations which include s-o coupling and 
also with VASP-GGA calculations with and without s-o interaction. The latter corroborate 
the general trend that s-o interaction decreases C.y by about 10%, a fact which corresponds 
to an increase in the TA phonon frequencies. The VASP calculations agree better with the 
experiment than the ABINIT ones but the excellent agreement with experimental data found 
for the VASP calculations without s-o interaction should, of course, be considered fortuitous. 
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The poor agreement of the ABINIT calculations with the experimental results corresponds 
to the poor agreement found in Fig. [6] for the TA-phonons along the F-X direction. The 
reason for the poor agreement found in this case is not known to us in detail. We presume 
it is due to inaccuracies in the ab initio pseudopotentials employed. 

At high temperatures contributions linear in temperature due to thermal expansion have 
been added to the ABINIT data (see insets in Fig. [13] and Fig. [T^ . While the linear term 
for HgTe (9.08x10"'^ J/molK^ T) is in good agreement with literature data, that for HgSe 
(~4x 10^^ J/mol T) is smaller by a factor of 3.5 than data reported in the literature.— Our 
experimental results are, however, in very good agreement with data reported by Gul'tyaev 
and coworkers.— 




FIG. 13: (color online) Cy/T'^ calculated for HgSe vs T with the VASP code (with (red dashed 
line) and without (blue solid line) ( s-o interaction) and (insert) with the ABINIT code (only 
with s-o interaction). The effect of the latter is similar to that found for the VASP calculations). 
Also, Cp/T^ measured by us and measurements by Gul'tyaev and Petrov and Khattak et al. are 
presented.—*^ The inset depicts Cp obtained by ABINIT with s-o interaction included. A 
contribution linear in temperature, 4x10"^ J/mol T, due to thermal expansion has been added 
to the calculated data. 

We would like to recall the fact that the s-o effects calculated here for materials with 
zincblende structure are opposite to those calculated in Romero et al.— for the lead chalco- 
genides, which crystallize in the rock salt structure. In order to investigate whether this 
difference is characteristic of the crystal structures, we calculated the dispersion relations of 
rock salt type HgS which, as shown in Fig. 9, is stable above 21.6 GPa (the experimental 
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FIG. 14: C^/T3 vs. T for HgTe calculated with both, ABINIT-LDA and VASP-GGA codes 
including and not including s-o interaction compared with our experimental data. The inset depicts 
C^,~ Cp obtained by ABINIT with s-o interaction included. A contribution linear in temperature, 



9.08x10-3 J/molK2 T (Ref. 
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), due to thermal expansion has been added to the calculated data. 



phase transition pressure is 20.5 GPai^). There are, of course, no measurements for HgS in 
the rock salt phase but we show in Fig. [15] the Ci,/T^ vs. T calculated with the VASP-GGA 
code. In order to stabilize the rock salt structure, we performed the calculations assuming a 
pressure of 25 GPa applied to the material. The s-o interaction induces a minute decrease 
of Cy/T^ (less than 1%) which can be assumed to be within the possible error of the cal- 
culations using different potentials. The change in the sign of the effect when going from 
zincblende to rock salt, found in the calculations going from the Hg chalcogenides to the 
Pb chalcogenides, is not confirmed by the calculations performed for HgS for the rock salt 
structure. 

The difference between zincblende and rocksalt, however, has the right sign. We thus 
conclude that the reason for the change in sign of the effect of s-o interaction, from the Hg 
to the Pb chalcogenides is likely to be related to the change in structure. 

The reader may wonder why the maximum of C^/T^ in Fig. [15] is smaller by a factor of 
seven than calculated for the zb phase. In order to ascertain the reason for this difference 
we calculated the phonon DOS for the rs phase. We found the TA peak to occur at 62 cm^^, 
a factor of 2.3 larger than that shown in Fig. [TWb) for the zb phase. This difference can 
account for the lower maximum in C.^,/T'^ calculated for the rs phase. 
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FIG. 15: (color online) Calculations, using the VASP-GGA code with (red) dashed line and without 
(solid line) s-o interaction, of HgS at a pressure of ^■25 GPa (rock salt structure). Notice that the 
s-o interaction lowers C^/T'^ by an amount that is insignificant with respect to uncertainties in the 
computational method. The inset has been added in order to clearly illustrate the small effect of 
s-o interaction. 

VII. CONCLUSIONS 

The original motivation of the present work was to perform accurate measurements of 
the specific heat of /3-HgS, HgSe and HgTe at low temperatures and to compare them with 
ah initio calculations of the band structures of these materials, thought to be semimetals 
with the /5-tin inverted valence band structure. A prerequisite was the calculation of the 
electronic band structures of these materials. While doing so, we found a number of thus 
far unknown details and their interpretation. We mention first that there is a significant 
hybridization of the d core electrons of Hg with the (chalcogen) pASke valence bands which 
results in a negative contribution to their s-o splitting at V. Consequently, for /3-HgS this s-o 
splitting even reverses its sign. Thus, the material, believed to be a semimetal, becomes a 
semiconductor with a gap similar to the magnitude of the F s-o splitting (~ 0.1 eV). The s-o 
splittings calculated for the three materials compare satisfactorily with experimental results, 
whereas the negative Eq gaps are larger in magnitude as a consequence of the so-called gap 
problem associated with LDA. The presence of s-o interaction in our calculations allows us 
to determine the coefficients of the linear terms of the Fg valence bands. They compare well 
with estimates based on a semiempirical interpolation formula used in the past for other 
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zincblende-type semiconductors. 

We next proceeded to the calculation of the phonon dispersion relations and comparison 
with experimental results based on Raman and INS spectroscopies. We performed these 
calculations using existing computer codes and investigated the effect of s-o interaction 
which turned out, in the case of the acoustic phonons, to have a sign opposite to that found 
for the rock-salt structure of the Pb chalcogenides. This effect resulted in changes by about 
10% in the values of the calculated maximum of C^/T^ vs. T: the s-o splitting lowers the 
value at this maximum, an effect which is opposite to that found for the rock-salt structure 
lead chalcogenides. An argument was found to correlate the sign of this effect with the 
structure of the materials by calculating C^/T^ for HgS in the rock salt structure. From the 
calculated phonon dispersion relations we obtained the phonon densities of states and their 
atomic projections. They were used for an analysis of measured phonon line widths and 
the corresponding effects of isotopic composition. We concluded that there must be a large 
contribution of the electron-phonon interaction across the Fg bands to the line widths of the 
Raman phonons of HgSe and HgTe. We also calculated two-phonon densities of states in 
order to interpret extant and future second order Raman spectra. Finally we investigated a 
few effects of hydrostatic pressure on the Raman phonons of these materials, in particular 
HgTe. We calculated the Griineisen parameters of LO and TO phonons and also of the 
transverse (Born) effective charge. The latter was found to be 7* = -0.65 for HgTe, a value 
similar to that found for most other zincblende type semiconductors, the sign of 7* has been 
interpreted as signaling a decrease of the ionicity with increasing pressure. 

Acknowledgments 

A.H.R. has been supported by CONACyT Mexico under project J-59853-F and by 
PROALMEX/DAAD. Further computer resources have been provided by CNS IPICYT 
Mexico. A. M. acknowledges the financial support from the Spanish MCYT under grants 
MAT2007-65990-C03-03, CSD2007-00045 and the supercomputer resources provides by the 
Red Espanola de Supercomputacion. We thank E. Briicher for the magnetic measurements 
and V. Duppel for the microprobe analysis of the natural /3-HgS crystals. We are also in- 
debted to K. Syassen for a critical reading of the manuscript and to N. E. Christensen for 



28 



pointing out to us references 



40 



44, 



45. 



* Corresponding author: E-mail r.kremer@fkf.mpg.de 



^ Cp represents the experimentally measured value of the heat capacity at constant pressure. The 

ah initio calculations yield i.e. the heat capacity at constant volume. In the temperature 

region considered here C„ w Cp. 
2 M. Sanati, S.K. Estreicher, and M.Cardona, Solid State Commun., 131, 229 (2004). 
^ A. Gibin, G. G. Devyatykh, A. V. Gusev, R. K. Kremer, M. Cardona, and H.-J. Pohl, Solid 

State Commun. 133, 569 (2005). 
^ R. K. Kremer, M. Cardona, E. Schmitt, J. Blumm, S.K. Estreicher, S. Sanati, M. Bockowski, 

I. Grzegory, T. Suski, and A. Jezowski, Phys. Rev. B 72, 075209 (2005). 
^ M. Cardona, R. K. Kremer, M. Sanati, S. K. Estreicher, T. R. Anthony, Solid State Commun. 

133, 465 2005. 

^ J. Serrano, R. K. Kremer, M. Cardona, G. Siegle, A. H. Romero, and R. Lauck, Phys. Rev. B 
73, 094303 (2006). 

^ A. H. Romero, M. Cardona, R. K. Kremer, R. Lauck, G. Siegle, J. Serrano, and X. C. Gonze, 

Phys. Rev. B 78, 224302 (2008). 
^ L. E. Dfaz-Sanchez, A. H. Romero, M. Cardona, R. K. Kremer, and X. Gonze, Phys. Rev. Lett. 

99, 165504 (2007). 

^ J. Serrano, R. K. Kremer, M. Cardona, G. Siegle, L. E. Diaz-Sanchez, and A. H. Romero, Phys. 

Rev. B 77, 054303 (2008). 
^° R. Shaltaf, X. Gonze, M. Cardona, R. K. Kremer, and G. Siegle, Phys. Rev. B 79, 075204 
(2009). 

S. H. Groves and W. Paul, Phys. Rev. Lett. 11, 194 (1963). 

R. Zallen and M. Slade, Solid State Commun. 8, 1291 (1970). Data for Aq and £"0 quoted in 
Tabled] were taken at 300 K. 

W. Szuszkiewicz, B. Witkowska, M. Jouanne and M. Balkanski, Materials Science Forum 182- 
184, 711 (1995). 

14 G. Masse, J. P. Aicardi and J. P. Leyris, J. Lumin. 17, 29 (1978). 
1^ S. Bloom and T. K. Berstresser, phys. stat. sol. 42, 191 (1970). 

29 



16 H. Overhof, phys. stat. sol., 43, 221 (1971). 

N. A. Cade and P. M. Lee, Solid State Commun. 56, 637 (1985). 
1^ M.Cardona, Phys. Rev. 129, 1068 (1963). 

1^ K. Shindo, A. Morita, and H. Kamimura, Proc. Phys. Soc. Japan 20, 2054 (1965). 

20 B. Tell and P. M. Bridenbaugh, Phys. Rev. B 12, 3330 (1975). 

21 G. Dresselhaus, Phys. Rev. 100, 580 (1955). 

22 M. Cardona, N. E. Christensen and G. Fasol, Phys. Rev. B 38, 1806 (1988). 

23 P. Hohenberg and W. Kohn, Phys. Rev. 136, B864 (1964). 

2^ X. Gonze, J.-M. Beuken, R. Caracas, F. Detraux, M. Fuchs, G.-M. Rignanese, L. Sindic, 
M. Verstraete, G. Zerah, F. Jollet, M. Torrent, A. Roy, M. Mikami, P. Ghosez, J.-Y. 
Raty, and D. C. Allan, Comput. Mater. Sci. 25, 478 2002 . ABINIT is a common project 
of the Universite Catholique de Louvain, Corning Incorporated, and other contributors 



http://www.pcpm.ucl.ac.be/ABINIT 



25 



G. Kresse and J. Furthmiiller, Comp. Mat. Sci. 6, 15 (1996), Phys. Rev. B 54, 11169 (1996) , 



more info of computer code VASP. See: http://cms.mpi.univie.ac.at/vasp 

26 W. Kohn and L. J. Sham, Phys. Rev. 140, A1133 (1965). 

27 D.M. Ceperly, B.J. Alder, Phys. Rev. Lett. 45, 566 (1980). 

28 J. P. Perdew, S. Burke, and M. Ernzerhof, Phys. Rev. Lett. 77, 3865 (1996). 

29 C. Hartwigsen, S. Goedecker, and J. Mutter, Phys. Rev. B 58, 3641 (1998). 
3° S. Goedecker, M. Teter, and J. Mutter, Phys. Rev. B 54, 1703 (1996). 

31 L. Kleinman, Phys. Rev. B 21, 2630 1980. 

32 L. A. Memstreet, C. Y. Fong, and J. S. Nelson, Phys. Rev. B 47, 4238 (1993). 

33 X. C. Gonze and C. Lee, Phys. Rev. B 55, 10355 1997. X. C. Gonze, Phys. Rev. B 55, 10337 
(1997). 

34 P. E. Blochl, Phys. Rev. B 50, 17953 (1994); G. Kresse and D. Joubert, Phys. Rev. B 59, 1758 
(1999). 

3^ F. J. Manjon, D. Errandonea, N. Garro, J. Pellicer-Porres, P. Rodriguez-Hernandez, S. Radescu, 
J. Lopez-Solano, A. Mujica, and A. Muhoz, Phys. Rev. B 74, 144111 (2006). 



36 K. Parlinski, computer code PMONON. See: http://wolf.ifj.edu.pl/phonon 
3^7 R. M. Pick, M. M. Cohen and R. M. Martin, Phys. Rev. B 1, 910 (1970). 



38 



F. Merman and S. Skillman, Atomic Energy Calculations, Prentice Hall, Englewood Cliffs, NJ, 



30 



1963. 

39 A. Delin, Phys. Rev. B65, 153205 (2002). 

^0 P. Carrier and S.-H. Wei, Phys. Rev. B 70, 035212 (2004). 

Landolt-Bornstein Tables, New Series Vol III/22a edited by O. Madelung (Springer 1986). 

E. S. Rittner and J. H. Schulman, J. Phys. Chem. 47, 537 (1943). 
^3 R. W. Godby, M. Schluter, and L. J. Sham Phys. Rev. Lett. 56, 2415 (1986); Phys. Rev. B 36, 

6497 (1987); Phys. Rev. B 37, 10159 (1988). 

44 M. Rohlfing and S. G. Louie, Phys. Rev. B57, R9392 (1998). 

45 A. Fleszar and W. Hanke, Phys. Rev. B71, 045207 (2005). 

46 M.P. Surh, M.F. Li, and S. G. Louie, Phys. Rev. B 43, 4286 (1991). 

4"^ Y. S. Kim, K. Hummer, and G. Kress, Phys. Rev. B 80, 035203 (2009). 
4^ W. Szuszkiewicz, B. Witkovska, M. Jouanne, Acta Phys. Polonica a 87, 415 (1995). 
49 W. Szuszkiewicz, K. Dybko, B. Hennion, M. Jouanne, C. Julien, E. Dynowska, J. Gbrecka, and 
B. Witkowska, J. Cryst. Growth 184-185, 1204 (1998). 

H. Kepa, T. Giebultowicz, B. Buras, B. Lebech, and K. Clausen, Phys. Scripta 25, 807 (1982); H. 
Kepa, H. Gebicki, T. Giebultowicz, B. Buras, B. Lebech, and K. Clausen, Solid State Commun. 
34, 211 (1980). 

51 A. Gobel, T. Ruf, J. H. Zhang, R. Lauck, and M. Cardona, Phys. Rev. B 59, 2749 (1999). 

52 M. P. Chamberlain, C. Trallero-Giner, and M. Cardona, Phys. Rev. B 51, 1680 (1995). 

53 P. G. Etchegoin, M. Cardona, R. Lauck, R. J. H. Clark, J. Serrano, A. R. Romero phys. stat. 
sol. (b), 245 1125 (2008). 

54 D. Olego and M. Cardona, Solid State Commun. 39, 1071 (1981). 

55 J. Lazewski, K. Parlinski, W. Szuszkiewicz, and B. Hennion, Phys. Rev B 67, 094305 (2003). 

56 A. Mooradian and T. C. Harman, in The Physics of Semimetals and Narrow Gap Semiconduc- 
tors, edited by D. L. Carter and R. T. Bate (Pergamon, Oxford, 1973) p. 297. 

5^ W. Szuszkiewicz, E. Dynowska, J. Gorecka, B. Witkowska, M. Jouanne, J. F. Morhange, C. 
Julien, and B. Hennion, phys. stat. sol. B 215, 93 (1999). 

58 J. F. Kielkopf, J.O.S.A 63, 987(1973). 

59 M. Cardona and M. L. W. Thewalt, Rev. Mod. Physics 77, 1173 (2005). 

60 A. Debernardi, Solid State Commun. 113, 1 (2000). 

61 M. Chandrasekhar, H. R. Chandrasekhar, M. Grimsditch, and M. Cardona, Phys. Rev. B 22, 

31 



4825 (1980). 

A. Mujica, A. Rubio, A. Munoz, and R. J. Needs, Rev. Mod. Physics 75, 863 (2003). 

R. J. Nelmes and M. I. McMahon, in Semiconductors and Semimetals edited by T. Suski and 

W. Paul (Academic, New York, 1998), Vol. 54, p. 145. 

A. San Miguel, N. G. Wright, M. I. McMahon, and R. J. Nelmes, Phys. Rev. B 51, 8731 (1995). 

E. Anastassakis and M. Cardona, in High Pressure in Semiconductor Physics II, edited by T. 

Suski and W. Paul (Academic, San Diego, 1998), p. 118. 

A. R. Gofii, K. Syassen, and M. Cardona, Phys. Rev. B 41, 10104 (1990). 

W. Schnelle and E. Gmelin, Thermochim. Acta, 391, 41 (2002). 

P. V. Gul'tyaev and A. V. Petrov, Sov. Phys. Solid State 1, 330 (1959). 

J. G. Collins, G. K. White, J. A. Birch, and T. F. Smith, J. Phys. C: Solid State Phys. 13, 1649 
(1980). 

G. D. Khattak, H. Akbarzadeh, and P. H. Keesom, Phys. Rev. B 23, 2911 (1981). 

H. Kobayashi, T. Nozue, T. Matsumura, T. Suzuki, and T. Kamimura, J. Phys.: Condens. 
Matter 11, 8673 (1999). 



32 



